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We extract parameters relevant for distinguishing among single-field inflation models from the 
Wilkinson Microwave Anisotropy Probe (WMAP) data set, and from a combination of the WMAP 
data and seven other Cosmic Microwave Background (CMB) experiments. We use only CMB data 
and perform a likelihood analysis over a grid of models including the full error covariance matrix. 
We find that a model with a scale-invariant scalar power spectrum (n = 1), no tensor contribution, 
and no running of the spectral index, is within the 1-er contours of both data sets. We then apply 
the Monte Carlo reconstruction technique to both data sets to generate an ensemble of inflationary 
potentials consistent with observations. None of the three basic classes of inflation models (small- 
field, large-field, and hybrid) are completely ruled out, although hybrid models are favored by the 
best-fit region. The reconstruction process indicates that a wide variety of smooth potentials for the 
inflaton are consistent with the data, implying that the first-year WMAP result is still too crude 
to constrain significantly either the height or the shape of the inflaton potential. In particular, the 
lack of evidence for tensor fluctuations makes it impossible to constrain the energy scale of inflation. 
Nonetheless, the data rule out a large portion of the available parameter space for inflation. For 
instance, we find that potentials of the form V = X(j) 4 are ruled out to 3a by the combined data set, 
but not by the WMAP data taken alone. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the fundamental ideas of modern cosmology 
is that there was an epoch early in the history of the 
universe when potential, or vacuum, energy dominated 
other forms of energy densities such as matter or radia- 
tion. During such a vacuum-dominated era the scale fac- 
tor grew exponentially (or nearly exponentially) in some 
small time. During this phase, dubbed inflation [jj, l^j , a 
small, smooth spatial region of size of order the Hubble 
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radius grew so large that it easily could encompass the 
comoving volume of the entire presently observable uni- 
verse. If the universe underwent such a period of rapid 
expansion, one can understand why the observed universe 
is homogeneous and isotropic to such high accuracy. 

One of the predictions of the simplest models of in- 
flation is a spatially flat Universe, i.e., Qtot — 1, with 
great precision. Inflation has also become the domi- 
nant paradigm for understanding the initial conditions 
for structure formation and for Cosmic Microwave Back- 
ground (CMB) anisotropy. In the inflationary picture, 
primordial density and gravity- wave (tensor) fluctuations 
are created from quantum fluctuations and "redshifted" 
out of the horizon during an early period of superlumi- 
nal expansion of the universe, where they are "frozen" 
as perturbations in the background metric [E El IE IE • 
Metric perturbations at the surface of last scattering are 
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observable as temperature anisotropics in the CMB. The 
first and most impressive confirmation of the inflation- 
ary paradigm came when the CMB anisotropics were 
detected by the Cosmic Background Explorer (COBE) 
satellite in 1992 [1 H [ljj. Subsequently, it became 
clear that the measurements of the spectrum of the CMB 
anisotropy can provide very detailed information about 
fundamental cosmological parameters 11] and other cru- 
cial parameters for particle physics. 

In the past few years, a number of balloon-borne and 
terrestrial experiments have mapped out the CMB an- 
gular anisotropies El El El EH El E3 El , revealing 
a remarkable agreement between the data and the infla- 
tionary predictions of a flat universe with a nearly scale- 
invariant spectrum of adiabatic primordial de nsity per- 
turbations, (see e.g., 

Despite the simplicity of the inflationary paradigm, the 
number of inflation models that have been proposed in 
the literature is enormous Q- This is true even if we 
limit ourselves to models with only one scalar field (the 
inflaton). With the previous data on CMB anisotropies 
from balloon and terrestrial experiments it has been pos- 
sible for the first time to place interesting constraints on 
the space of possible inflation models [28l 12^. l30| . How- 
ever, the quality of the data were not good enough to 
rule out entire classes of models. A boost along these 
lines has been very recently provided by the data of the 
Wilkinson Microwave Anisotropy Probe (WMAP) mis- 
sion, which has marked the beginning of the precision 
era of the CMB measurements in space 0, |3^. The 
WMAP collaboration has produced a full-sky map of the 
angular variations in the microwave flux, in particular the 
cosmic microwave background, with unprecedented accu- 
racy. WMAP data support the inflationary mechanism 
for the generation of curvature superhorizon fluctuations 
and provide a strong bound on the possible admixture 
of isocurvature modes [H- Furthermore, consistent with 
the simplest single- field models of inflation [34|], no evi- 
dence of nongaussianity is found [35| . 

The goal of this paper is to use the WMAP data to 
discriminate among the various single-field inflationary 
models. To obtain some indication of the robustness 
of our analysis, we also consider a data set consisting 
of WMAP augmented with several other CMB experi- 
ments. For single-field inflation models, the relevant pa- 
rameter space for distinguishing among models is defined 
by the scalar spectral index n, the ratio of tensor to scalar 
fluctuations r, and the running of the scalar spectral in- 
dex dn/dhik. We employ Monte Carlo reconstruction, 
a stochastic method for "inverting" observational con- 
straints to determine an ensemble of inflationary poten- 
tials compatible with observation |3rl W\ . In addition to 
encompassing a broader set of models than usually con- 
sidered (large-field, small-field, hybrid and linear mod- 
els) , Monte Carlo reconstruction makes it possible easily 
to incorporate constraints on the running of the spectral 
index as well as to include effects to higher order in slow 
roll. 



Since studies on the implications of WMAP data for 
inflation 33, 38] have already appeared, we briefly men- 
tion the different elements between our analysis and 
others. (We will elaborate on these differences later.) 
The WMAP collaboration analysi s Ipl ] included WMAP 
data, additional CMB data (CBI [l|| and ACBAR [H), 
large-scale structure data (2dFGRS H3), as well as 
Lyman-a power spectrum data [41| . They used a Markov 
Chain Monte Carlo technique to explore the likelihood 
surface. Barger et al. |H, considered WMAP data only, 
but with a top-hat prior on the Hubble constant h (Hp = 
100 h km sec" 1 Mpc" 1 ) from the HST key project [42] . 
Also, Barger et al. did not consider a running of the scalar 
spectral index. We only consider CMB data. We first 
analyze just the WMAP results. We then analyze the 
WMAP data set in conjunction with other CMB data 
sets (BOOMERanG-98 MAXIMA-1 0, DASI Jll, 
CBI E3, ACBAR H, VSAE gj|, and Archeops [lTf). 
We employ a grid of models in the likelihood analysis, 
which differs from the method used by the WMAP team. 

The paper is organized as follows. In Sec. II we dis- 
cuss single-field inflation and the relevant observables in 
more detail. In Sec. Ill we discuss the inflationary model 
space, and in Sec. IV we describe the Monte Carlo re- 
construction technique. Section V describes the methods 
used for the CMB analysis. In Sec. VI we present the 
constraints from the CMB anisotropy data sets. In Sec. 
VII we present our conclusions. 



II. SINGLE-FIELD INFLATION AND THE 
INFLATIONARY OBSERVABLES 

In this section we briefly review scalar field models 
of inflationary cosmology, and explain how we relate 
model parameters to observable quantities. Inflation, in 
its most general sense, can be defined to be a period 
of accelerating cosmological expansion during which the 
universe evolves toward homogeneity and flatness. This 
acceleration is typically a result of the universe being 
dominated by vacuum energy, with an equation of state 
p ~ —p. Within this broad framework, many specific 
models for inflation have been proposed. We limit our- 
selves here to models with "normal" gravity (i.e., general 
relativity) and a single order parameter for the vacuum, 
described by a slowly rolling scalar field 4>, the infla- 
ton. These assumptions are not overly restrictive; the 
most widely studied inflation models fall within this cat- 
egory, including Linde's "chaotic" inflation scenario |4tS| . 
inflation from pseudo Nambu-Goldstone bosons ("nat- 
ural" inflation 47]), dilaton-like models involving ex- 
ponential potentials (power-law inflation), hybrid infla- 
tion [481 l49l l50j . and so forth. Other models, such as 
Starobinsky's R 2 model [&l| an d versions of extended in- 
flation, can, through a suitable transformation, be viewed 
in terms of equivalent single- field models. Of course in 
single-field models of inflation, the inflaton "field" need 
not be a fundamental field at all. Also, some "single- 
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field" models require auxiliary fields. Hybrid inflation 
models 0, |5(j, for example, require a second field 
to end inflation. What is significant is that the infla- 
tionary epoch be described by a single dynamical order 
parameter, the inflaton field. 

A scalar field in a cosmological background evolves 
with an equation of motion 



4> + 3H<j)+V' ((f)) =0. 



(1) 



The evolution of the scale factor is given by the scalar 
field dominated FRW equation, 
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GeV is the Planck mass, 
and we have assumed a flat Friedmann- Robertson- Walker 
metric, 



= diag(l, 



a 2 ), 



(3) 



where a 2 (t) is the scale factor of the universe. Inflation is 
defined to be a period of accelerated expansion, a > 0. A 
powerful way of describing the dynamics of a scalar field- 
dominated cosmology is to express the Hubble parameter 
as a function of the field <f>, H = H(<f>), which is consistent 
provided <j> is monotonic in time. The equations of motion 
become [52. 1531 l54l l55| : 
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These are completely equivalent to the second-order 
equation of motion in Eq. QJ. The second of the above 
equations is referred to as the Hamilton- Jacobi equation, 
and can be written in the useful form 



H 2 (4>) 



/ 8vr 



where e is defined to be 
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The physical meaning of e can be seen by expressing Eq. 
© as 



= H 2 ^) [l-e(0], 



(7) 



so that the condition for inflation (a/a) > is given by 
e < 1. The scale factor is given by 
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where the number of e-folds N is 
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To create the observed flatness and homogeneity of the 
universe, we require many e-folds of inflation, typically 
N ~ 60. This figure varies somewhat with the details 
of the model. We can relate a comoving scale k in the 
universe today to the number of e-folds N before the end 
of inflation by [Hij 



N(k) = 62-ln — - — In 
a H 
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(10) 

Here Vk is the potential when the mode leaves the hori- 
zon, V e is the potential at the end of inflation, and prh 
is the energy density after reheating. Scales of order the 
current horizon size exited the horizon at (k) ~ 60. 
Since this number depends, for example, on the details 
of reheating, we will allow N to vary within the range 
40 < N < 70 for any given model in order to consider 
the most general case. (Dodelson and Hui have recently 
argued that the value of N corresponding to the current 
horizon size can be no larger than 60 |57| , so in this sense 
we are being more general than is necessary.) 

We will frequently work within the context of the slow 
roll approximation |58l [EjJ , which is the assumption that 
the evolution of the field is dominated by drag from the 
cosmological expansion, so that (j> ~ and 



3iT 



(11) 



The equation of state of the scalar field is dominated by 
the potential, so that p ~ — p, and the expansion rate is 
approximately 



H 



I 8tt 
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(12) 



The slow roll approximation is consistent if both the slope 
and curvature of the potential are small, V', V" <C V. 
In this case the parameter e can be expressed in terms of 
the potential as 



'pi 



(13) 
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We will also define a second "slow roll parameter" rj by: 
'H» (0)' 
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Slow roll is then a consistent approximation for e, r/ <C 1. 
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Inflation models not only explain the large-scale ho- 
mogeneity of the universe, but also provide a mecha- 
nism for explaining the observed level of inhomogene- 
ity as well. During inflation, quantum fluctuations on 
small scales are quickly redshifted to scales much larger 
than the horizon size, where they are "frozen" as per- 
turbations in the background metric 

@,aill0 The 

metric perturbations created during inflation are of two 
types: scalar, or curvature perturbations, which couple to 
the stress-energy of matter in the universe and form the 
"seeds" for structure formation, and tensor, or gravita- 
tional wave perturbations, which do not couple to matter. 
Both scalar and tensor perturbations contribute to CMB 
anisotropy. Scalar fluctuations can also be interpreted as 
fluctuations in the density of the matter in the universe. 
Scalar fluctuations can be quantitatively characterized 
by perturbations P-n in the intrinsic curvature scalar. As 
long as the equation of state e is slowly varying, 1 the cur- 
vature perturbation can be shown to be |6]J, |62j, l63l l64| 



The ratio of tensor to scalar modes is then 



,-,1/2 
R 



Z7T 
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(15) 
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The fluctuation power spectrum is in general a function 
of wavenumber k, and is evaluated when a given mode 
crosses outside the horizon during inflation, k = aH. 
Outside the horizon, modes do not evolve, so the ampli- 
tude of the mode when it crosses back inside the horizon 
during a later radiation- or matter-dominated epoch is 
just its value when it left the horizon during inflation. 
The spectral index n for P-r is defined by 



n-l = 



dlnPn 
dhik 



(16) 



so that a scale-invariant spectrum, in which modes have 
constant amplitude at horizon crossing, is characterized 
by n = 1 . Some inflation models predict run ning of the 
spectral index with scale [HH Ea EE E3, IZljT^, E3l 
or even sharp features in the power spectrum |74| . We 
will consider the running of the spectral index in more 
detail on Sec. IIVI 

Instead of specifying the fluctuation amplitude directly 
as a function of k, it is often convenient to specify it as 
a function of the number of e-folds N before the end 
of inflation at which a mode crossed outside the hori- 
zon. Scales of interest for current measurements of CMB 
anisotropy crossed outside the horizon at N ~ 60, so that 
P-jz is conventionally evaluated at P-r (N ~ 60). 

The power spectrum of tensor fluctuation modes is 
given by [UElESElEl 
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1 This assumption is not identical to the assumption of slow roll 
(see, e.g., Ref. l60l ). although in most cases it is equivalent. 
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so that tensor modes are negligible for e <C 1. Tensor 
and scalar modes both contribute to CMB temperature 
anisotropy. If the contribution of tensor modes to the 
CMB anisotropy can be neglected, normalization to the 



COBE four-year data gives [80J, |81| P 



,1/2 
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4.8 x 10- 



In the next section, we will describe the predictions of 
various models in this parameter space. 



III. THE INFLATIONARY MODEL SPACE 

To summarize the results of the previous section, infla- 
tion generates scalar (density) and tensor (gravity wave) 
fluctuations which are generally well approximated by 
power laws: 



Pn (k) cx k r 



P T (k) cx k r ' 



(19) 



In the limit of slow roll, the spectral indices n and tit 
vary slowly or not at all with scale. We can write the 
spectral indices n and nx to lowest order in terms of the 
slow roll parameters e and r\ as |64| : 



n 
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(20) 



The tensor spectral index is not an independent param- 
eter, but is proportional to the tensor/scalar ratio, given 
to lowest order in slow roll by 



Pt 

' 2e = - 2 P^- 



(21) 



This is known as the consistency relation for inflation. 
(This relation holds only for single-field inflation, and 
weakens to an inequality for inflation involving multiple 
degrees of freedom j82,|83|,|84(.) A given inflation model 
can therefore be described to lowest order in slow roll 
by three independent parameters, Pr, Pt, and n. If 
we wish to include higher-order effects, we have a fourth 
parameter describing the running of the scalar spectral 
index, dn/d In k. 

The tensor/scalar ratio is frequently expressed as a ra- 
tio of their contributions to the CMB quadrupole, 



^Tensor 
^Scalar 



(22) 



The relation between r and ratio of amplitudes in the 
primordial power spectra Pt/Pk depends on the back- 
ground cosmology, in particular the densities of matter 
(O m ) and cosmological constant (Oa)- For the currently 
favored values of f2 m ~ 0.3 and J7a — 0.7, the relation is 
approximately 



r ~ lOe, 



(23) 
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to lowest order in slow roll. Conventions for the normal- 
ization of this parameter vary widely in the literature. In 
particular, Peiris et al. [23] use r ~ 16e. 

Calculating the CMB fluctuations from a particular 
inflationary model reduces to the following basic steps: 
(1) from the potential, calculate e and r\. (2) From e, 
calculate ./V as a function of the field <fi. (3) Invert N (</>) 
to find 4>m- (4) Calculate Pji, n, and Pt as functions of 
(f>, and evaluate them at <fi — <fipf. For the remainder of 
the paper, all parameters are assumed to be evaluated at 

4> = 4>n- 

Even restricting ourselves to a simple single-field infla- 
tion scenario,the number of models available to choose 
from is large [2] • It is convenient to define a general clas- 
sification scheme, or "zoology" for models of inflation. 
We divide models into three general types: large-field, 
small-field, and hybrid, with a fourth classification, lin- 
ear models, serving as a boundary between large- and 
small-field. A generic single-field potential can be char- 
acterized by two independent mass scales: a "height" A 4 , 
corresponding to the vacuum energy density during in- 
flation, and a "width" fi, corresponding to the change in 
the field value Aip during inflation: 



V((f>)=A 4 f 



(24) 



Different models have different forms for the function /. 
The height A is fixed by normalization, so the only free 
parameter is the width /i. 

With the normalization fixed, the relevant parameter 
space for distinguishing between inflation models to low- 
est order in slow roll is then the r — n plane. (To next or- 
der in slow-roll parameters, one must introduce the run- 
ning of n.) Different classes of models are distinguished 
by the value of the second derivative of the potential, or, 
equivalently, by the relationship between the values of 
the slow-roll parameters e and r\. 2 Each class of mod- 
els has a different relationship between r and n. For a 
more detailed discussion of these relations, the reader is 
referred to Refs. [H|. 

First order in e and r\ is sufficiently accurate for the 
purposes of this Section, and for the remainder of this 
Section we will only work to first order. The generaliza- 
tion to higher order in slow roll will be discussed in Sec. 

El 



A. Large-field models: — e < r\ < e 

Large-field models have inflaton potentials typical of 
"chaotic" inflation scenarios |4g , in which the scalar field 



2 The designations "small-field" and "large-field" can sometimes 
be misleading. For in stan ce, both the B? model l51l and the 
"dual inflation" model 85] are characterized by A(f> ~ Tipi, but 
are "small-field" in the sense that r\ < < e, with n < 1 and 
negligible tensor modes. 



is displaced from the minimum of the potential by an 
amount usually of order the Planck mass. Such models 
are characterized by V" ((f)) > 0, and — e < r) < e. The 
generic large-field potentials we consider are polynomial 
potentials V ((f)) = A 4 ((f>/n) p , and exponential potentials, 
V ((f)) — A 4 exp (4>/n). For the case of an exponential 
potential, V ((f)) cx exp ((f)/ fi), the tensor/scalar ratio r is 
simply related to the spectral index as 



r = 5 (1 — n) . 



(25) 



This result is often incorrectly generalized to all slow- 
roll models, but is in fact characteristic only of power- 
law inflation. For inflation with a polynomial potential, 
V ((f>) oc 4> p , we again have r cx 1 — n, 



r = 5 



P 



P- 



(1 



(26) 



so that tensor modes are large for significantly tilted spec- 
tra. We will be particularly interested in models with 
p = 4 as a test case for our ability to rule out models. 
For p = 4, the observables arc given in terms of the num- 
ber of e-folds N by 
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l-n = 



N + 1' 
3 

N + l' 



(27) 



B. Small-field models: r] < — e 



Small-field models are the type of potentials that arise 
naturally from spontaneous symmetry breaking (such as 
the original models of "new" inflation [58|,|59|) and from 
pseudo Nambu-Goldstone modes (natural inflation [4jj). 
The field starts from near an unstable equilibrium (taken 
to be at the origin) and rolls down the potential to a 
stable minimum. Small-field models are characterized by 
V" ((f)) < and n < —e. Typically e (and hence the 
tensor amplitude) is close to zero in small-field models. 
The generic small-field potentials we consider are of the 
form V (4>) = A 4 [1 — (<f>/ fi) p ], which can be viewed as a 
lowest-order Taylor expansion of an arbitrary potential 
about the origin. The cases p = 2 and p > 2 have very 
different behavior. For p = 2, 



r = 5(1 - n) exp [-1 — JV (1 — ri)] 



(28) 



where N is the number of e-folds of inflation. For p > 2, 
the scalar spectral index is 



i 2 fp- 1 
n ~ 1 — — 



N \p- 2 J ' 



(29) 



independent of r. Assuming /i < mp\ results in an upper 
bound on r of 



r < 5- 
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(30) 
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C. Hybrid models: < e < n 



IV. MONTE CARLO RECONSTRUCTION 



The hybrid scenario |48|, |49|, |50| frequently appears in 
models which incorporate inflation into supersymmctry. 
In a typical hybrid inflation model, the scalar field re- 
sponsible for inflation evolves toward a minimum with 
nonzero vacuum energy. The end of inflation arises as a 
result of instability in a second field. Such models are 
characterized by V" (<f) > and < e < rj. We con- 
sider generic potentials for hybrid inflation of the form 
V (0) = A 4 [1 + 0/Ai) p ] • The field value at the end of 
inflation is determined by some other physics, so there is 
a second free parameter characterizing the models. Be- 
cause of this extra freedom, hybrid models fill a broad 
region in the r — n plane (see Fig.^l. There is, however, 
no overlap in the r — n plane between hybrid inflation 
and other models. The distinguishing feature of many 
hybrid models is a blue scalar spectral index, n > 1. This 
corresponds to the case n > 2e. Hybrid models can also 
in principle have a red spectrum, n < 1. 



D. Linear models: 77 



In this section we describe Monte Carlo reconstruc- 
tion, a stochastic method for "inverting" observational 
constraints to determine an ensemble of inflationary po- 
tentials compatible with observation. The method is de- 
scribed in more detail in Refs. [3(1 l37j . In addition to 
encompassing a broader set of models than we consid- 
ered in Sec. IIIII Monte Carlo reconstruction allows us 
easily to incorporate constraints on the running of the 
spectral index dn/dh\k as well as to include effects to 
higher order in slow roll. 

We have defined the slow roll parameters e and i] in 
terms of the Hubble parameter H ((f)) as 



47T 
2 
PI 

4"7T 
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H{4>) 
H"{<P) 
H(cb) 



(32) 



These parameters are simply related to observables r ~ 
lOe, and n—1 ~ 4e— 2-q to first order in slow roll. (We dis- 
cuss higher order expressions for the observables below.) 
Taking higher derivatives of H with respect to the field, 
we can define an infinite hierarchy of slow roll parameters 



Linear models, V ((f)) oc 0, live on the boundary be- 
tween large-field and small-field models, with V" ((f)) = 
and n = — e. The spectral index and tensor/scalar ratio 
are related as: 



(1-n). 



(31) 



This enumeration of models is certainly not exhaus- 
tive. There are a number of single-field models that do 
not fit well into this scheme, for example logarithmic 
potentials V ((f)) cx ln(^) typical of supersymmetry 0. 
Another example is potentials with negative powers of 
the scalar field V ((f)) oc (f>~ p used in intermediate infla- 
tion |86| and dynamical supersymmetric inflation |69Ll7l| . 
Both of these cases require an auxiliary field to end infla- 
tion and are more properly categorized as hybrid models, 
but fall into the small-field region of the r — n plane. 
However, the three classes categorized by the relation- 
ship between the slow-roll parameters as — e < 77 < e 
(large-field), n < — e (small-field, linear), and < e < 77 
(hybrid), cover the entire r — n plane and are in that 
sense complete. 3 Figure 1 J2§| shows the r — n plane 
divided into regions representing the large field, small- 
field and hybrid cases. Figure 2 shows a "zoo plot" of 
the particular potentials considered here plotted on the 
r — n plane. 



3 Ref. 29] incorrectly specified < r\ < e for large- field and r\ < 
for small-field. 
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Here we have chosen the parameter a = 2rj — 4e ~ n 
to make comparison with observation convenient. 

It is convenient to use N as the measure of time during 
inflation. As above, we take t e and 4> e to be the time and 
field value at end of inflation. Therefore, N is defined 
as the number of e-folds before the end of inflation, and 
increases as one goes backward in time (dt > =>■ dN < 
0): 

d d 
dN 



d\na 2y / 7r 



mpi r d 
ve-77, 



(34) 



where we have chosen the sign convention that \[i has 
the same sign as H' ((f)) : 



TOpl H' 



(35) 



20? H 

Then e itself can be expressed in terms of H and N simply 

as, 

1 dH 

HdN 

Similarly, the evolution of the higher order parameters 
during inflation is determined by a set of "flow" equations 



e. 



(36) 



de 

In 



e (a + 2e) . 
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da 
dN 

dN 



-5ecr - 12e 2 + 2 ( 2 A H ) , 
£- 1 



2)e 



(%)+ £+1 A H .(37) 



The derivative of a slow roll parameter at a given order 
is higher order in slow roll. A boundary condition can be 
specified at any point in the inflationary evolution by se- 
lecting a set of parameters e, cr, 2 A H , ... for a given value 
of TV. This is sufficient to specify a "path" in the in- 
flationary parameter space that specifies the background 
evolution of the spacetime. Taken to infinite order, this 
set of equations completely specifies the cosmological evo- 
lution, up to the normalization of the Hubble parameter 
H. Furthermore, such a specification is exact, with no 
assumption of slow roll necessary. In practice, we must 
truncate the expansion at finite order by assuming that 
the 1 Ah are all zero above some fixed value of I. We 
choose initial values for the parameters at random from 
the following ranges: 

N = [40,70] 

e = [0,0.8] 

cr = [-0.5,0.5] 

2 Ah = [-0.05,0.05] 

H = [-0.025,0.025], 



3 A 



M+1 A H = 0. 



(38) 

Here the expansion is truncated to order M by setting 
+1 Ah = 0. In this case, we still generate an exact solu- 
tion of the background equations, albeit one chosen from 
a subset of the complete space of models. This is equiv- 
alent to placing constraints on the form of the potential 
V (4>) , but the constraints can be made arbitrarily weak 
by evaluating the expansion to higher order. For the pur- 
poses of this analysis, we choose M — 5. The results are 
not sensitive to either the choice of order M (as long as 
it is large enough) or to the specific ranges from which 
the initial parameters are chosen. 

Once we obtain a solution to the flow equations 
[e(A^), a(N), £ Xn(N)], we can calculate the predicted val- 
ues of the tensor/scalar ratio r, the spectral index n, and 
the "running" of the spectral index dn/d\nk. To low- 
est order, the relationship between the slow roll param- 
eters and the observables is especially simple: r = 10c, 
n — 1 = cr, and dn/dhik = 0. To second order in slow 
roll, the observables are given by 0, |§3 > 

r = lOe [1 - C (cr + 2e)] , (39) 

for the tensor/scalar ratio, and 

n - 1 = cr - (5 - 3C) e 2 - J (3 - 5C) <re+ i (3 - C) ( 2 A H ) 

(40) 

for the spectral index. The constant C = 4(ln2+7) — 5 = 
0.0814514, where 7 ~ 0.577 is Euler's constant. 4 Deriva- 



tives with respect to wavenumber k can be expressed in 
terms of derivatives with respect to N as |9(J 



d 

In 



(i 



d\nk 



(41) 



The scale dependence of n is then given by the simple 
expression 



dn 
c?ln k 



1 



dn 
dN' 



(42) 



which can be evaluated by using Eq. an d the flow 



equations. For example, for the case of V 
observables to lowest order are 

10 



the 



N + 1 



dn 
dhik 



N - 



- 1 
3 



N(N + 1) 



(43) 



The final result following the evaluation of a particular 
path in the M-dimensional "slow roll space" is a point in 
"observable parameter space," i.e., (r, n, dn/dhik), cor- 
responding to the observational prediction for that par- 
ticular model. This process can be repeated for a large 
number of models, and used to study the attractor be- 
havior of the inflationary dynamics. In fact, the models 
cluster strongly in the observable parameter space |3fij |. 
Figure shows an ensemble of models generated stochas- 
tically on the (r, n) plane, along with the predictions of 
the specific models considered in Sec. IIIII 

Figure 0] shows an ensemble of models generated 
stochastically on the (n, dn/d In k) plane. As one can 
see, and contrary to what commonly believed, there are 
single-field models of inflation which predict a significant 
running of the spectral index. The same can be appre- 
ciated in Fig. El where we plot an ensemble of models 
generated stochastically on the (r, dn/dhik) plane. 

The reconstruction method works as follows: 

1. Specify a "window" of parameter space: e.g., cen- 
tral values for n — 1, r, or dn/dhik and their asso- 
ciated error bars. 

2. Select a random point in slow roll space, [e, 77, 1 Xb], 
truncated at order M in the slow roll expansion. 

3. Evolve forward in time (dN < 0) until either (a) 
inflation ends (e > 1), or (b) the evolution reaches 
a late-time fixed point (e = ^Ah = 0, cr = const.). 

4. If the evolution reaches a late-time fixed point, cal- 
culate the observables r, n— 1, and dn/dhik at this 
point. 



4 Some earlier papers I '>(ll S7l . due to a long unnoticed typographic 



error in Ref. l87l . used an incorrect value for the constant C, 
given by C = 4(m2 + 7) = 5.0184514. The effect of this error is 
significant at second order in slow roll. 
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5. If inflation ends, evaluate the flow equations back- 
ward N e-folds from the end of inflation. Calculate 
the observable parameters at that point. 

6. If the observable parameters lie within the speci- 
fied window of parameter space, compute the po- 
tential and add this model to the ensemble of "re- 
constructed" potentials. 

7. Repeat steps 2 through 6 until the desired number 
of models have been found. 

The condition for the end of inflation is that e = 1. 
Integrating the flow equations forward in time will yield 
two possible outcomes. One possibility is that the con- 
dition e = 1 may be satisfied for some finite value of N, 
which defines the end of inflation. We identify this point 
as N — so that the primordial fluctuations are actu- 
ally generated when N ~ 60. Alternatively, the solution 
can evolve toward an inflationary attractor with r = 
and n > 1, in which case inflation never stops. 5 In re- 
ality, inflation must stop at some point, presumably via 
some sort of instability, such as the "hybrid" inflation 
mechanism 0, 0, [5(| . Here we make the simplifying 
assumption that the observables for such models are the 
values at the late-time attractor. 

Given a path in the slow roll parameter space, the form 
of the potential is fixed, up to normalization |37j, ISll 1921 
l93| . The starting point is the Hamilton- Jacobi equation, 



V(<f>) 



/3m 2 



pi 



V 87 



l-ge(0 



(44) 



We have e(N) trivially from the flow equations. In order 
to calculate the potential, we need to determine H(N) 
and <p(N). With e known, H(N) can be determined by 
inverting the definition of e, Eq. 136fl . Similarly, </>(N) 
follows from the first Hamilton- Jacobi equation Q : 



m PL 



dN 2^ 



(45) 



Using these equations and Eq. (|44|l . the form of the po- 
tential can then be fully reconstructed from the numerical 
solution for e(N) . The only necessary observational input 
is the normalization of the Hubble parameter H, which 
enters the above equations as an integration constant. 
Here we use the simple condition that the density fluctu- 
ation amplitude (as determined by a first-order slow roll 
expression) be of order 10 -5 , 



Sp^l H 1 
p 2ir mpi yfe 
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-5 



(46) 



A more sophisticated treatment would perform a full nor- 
malization to the COBE CMB data [Hill. The value 01 
the field, 0, also contains an arbitrary, additive constant. 



See Ref. Litl for a detailed discussion of the fixed-point structure 
of the slow roll space. 



V. CMB ANALYSIS 

Our analysis method is based on the computation 
of a likelihood distribution over a fixed grid of pre- 
computed theoretical models. We restrict our analy- 
sis to a flat, adiabatic, A-CDM model template com- 
puted with CMBFAST ([9(|), sampling the parameters 
as follows: fl c dmh 2 = uj c dm = 0.01, ...0.25, in steps of 
0.01; n b h 2 = u b = 0.009, ...,0.028, in steps of 0.001 and 
tt A = 0.5, ...,0.95, in steps of 0.05. The value of the 
Hubble constant is not an independent parameter, since: 



h = 



Ucdm + <^b 



(47) 



and we use the further prior: h = 0.72 ± 0.15. We allow 
for a reionization of the intergalactic medium by vary- 
ing the Compton optical depth parameter r c in the range 
t c = 0.05, 0.30 in steps of 0.05. Our choice of the above 
parameters is motivated by Big Bang Nucleosynthesis 
bounds on ujh (both from D |93 and 4 He+ 7 Li 98]), from 
supernovae (|99j|) and galaxy clustering observations (see 
e.g., |l00jp . and by the WMAP temperature-polarization 
cross-correlation data, which indicate an optical depth 
r = 0.17 ± 0.04 (|23). Our choice for an upper limit of 
t c < 0.30 is very conservative respect to the max imum 
values expected in numerical simulations (see e.g., [I0lj) 
even in the case of non-standard reionization processes. 
From the grid above, we only consider models with an 
age of the universe in excess of 11 Gyr. Variations in 
the inflationary parameters n, r and dn/dhik are not 
computationally relevant, and for the range of values we 
considered they can be assumed as free parameters. The 
tensor spectral index n t is determined by the consistency 
relation. 

For the WMAP data we use the recent temperature 
and cross polarization results from Ref. |3l| and compute 
the likelihood £ WM AP for each theoretical model as ex- 
plained in Ref. 102j, using the publicly available code on 
the LAMBDA web site (http://lambda.gsfc.nasa.gov/). 

We further include the results from seven other ex- 
periments: BOOMERanG-98 E^K MAXIMA-1 0, 
DASI HILCBI [3, ACBAR [33, VSAE 0, and 
Archeops |17| . The expected theoretical Gaussian signal 
inside the bin C^ 1 is computed by using the publicly avail- 
able window functions and lognormal prefactors as in Ref. 
[l03] . The likelihood C pre ~ WMAP from this dataset and 
for a given theoretical model is defined by 



21n£ pr 



e-WMAP 



{C^-C% X )M BB ,{C%^C% X , 



(48) 

where Mbb 1 is the Gaussian curvature of the likelihood 
matrix at the peak and C e B x is the experimental signal 
in the bin. We consider 7%, 10%, 4%, 5%, 5%, 5% 
and 5% Gaussian distributed calibration errors for the 
Archeops, BOOMERanG-98, DASI, MAXIMA-1, VSAE, 
ACBAR, and CBI, experiments respectively and include 
the beam uncertainties u sing t he analytical marginaliza- 
tion method presented in |l04j | . We use as combined like- 
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lihood just the normalized product of the two likelihood 
distributions: C ~ C WMAP x ff^™" 

In order to constrain a set of parameters x we marginal- 
ize over the values of the remaining "nuisance" parame- 
ters y. This yields the marginalized likelihood distribu- 
tion 

C(x) ee P{x\C B ) = J C(x, y)dy. (49) 

In the next section we will present constraints in sev- 
eral two-dimensional planes. To construct plots in two 
dimensions, we project (not marginalize) over the third, 
"nuisance" parameter. For example, likelihoods in the 
(r, n) plane are calculated for given choice of r and n by 
using the value of dn/dlnk which maximizes the like- 
lihood function. The error contours are then plotted 
relative to likelihood falloffs of 0.17, 0.018 and 0.0035 
as appropriate for la, 2a, and 3<r contours of a three- 
dimensional Gaussian. In effect we are taking the shadow 
of the three-dimensional error contours rather than a 
slice, which makes clear the relationship between the er- 
ror contours and the points generated by Monte Carlo. 
All likelihoods used to constrain models are calculated 
relative the the full three-dimensional likelihood function. 



VI. RESULTS 

We will plot the likelihood contours obtained from our 
analysis on three different planes: dn/dlnk vs. n, r vs. 
n and r vs. dn/dlnk. Presenting our results on these 
planes is useful for understanding the effects of theoreti- 
cal assumptions and/or external priors. 

We do this in Fig. for two cases: the WMAP 
dataset alone (left column), and WMAP plus the addi- 
tional CMB experiments BOOMERanG-98, MAXIMA- 
1, DASI, CBI, ACBAR, VSAE, and Archeops (right col- 
umn). By analyzing these different datasets we can check 
the consistency of the previous experiments with WMAP. 
The dots superimposed on the likelihood contours show 
the models sampled by the Monte Carlo reconstruction. 

In the top row of Fig. H3 we show the 68%, 95%, and 
99% likelihood contours on the n vs. dn/dlnk plane (refer 
to the end of Sec. for a discussion of the method used 
to plot the contours.) The pivot scale k$ is feo = 0.002/i 
Mpc" 1 . As we can see, both datasets are consistent with 
a scale invariant n — 1 power law spectrum with no fur- 
ther scale dependence (dn/dlnk = 0). A degeneracy is 
also evident: an increase in the spectral index n is equiv- 
alent to a negative scale dependence (dn/dlnk < 0). We 



6 We do not take into account correlations in the variance of the 
experimental data introduced by observations of the same por- 
tions of the sky (such as in the case of WMAP and Archeops, for 
example) . We have verified that such correlations have negligible 
effect by removing one experiment at a time and testing for the 
stability of our results. 



emphasize, however, that this beahavior depends strictly 
on the position of the pivot scale k : choosing k = 0.05h 
Mpc -1 would change the direction of the degeneracy. 
Models with n ~ 1.1 need a negative running at about the 
3a level. It is interesting also to note that models with 
lower spectral index, n ~ 0.9, are in better agreement 
with the data with a zero or positive running. For n < 1, 
the running is bounded by 0.005 > dn/dlnk > —0.025 
at la. 

In the center row of Fig. [fj] we plot the 68%, 95%, and 
99% likelihood contours on the r vs. n plane. As we can 
see, the present data only weakly constrain the presence 
of tensor modes, although a gravity wave component is 
not preferred. Models with n < 0.9 must have a neg- 
ligible tensor component, while models with n > 1 can 
have r larger than 0.4 (2a C.L.). However, as we can 
see from the bottom row of Fig.El there is no correlation 
between the tensor component and the running of the 
scalar index. 

Figures 03 and show a subset of 300 reconstructed 
potentials selected from the sampled set. Note in par- 
ticular the wide range of inflationary energy scales com- 
patible with the observational constraint. This is to be 
expected, since there was no detection of tensor modes in 
the WMAP data, which would be seen here as a detection 
of a nonzero tensor/scalar ratio r. In addition, the shape 
of the inflationary potential is also not well constrained 
by WMAP. 

Figure El shows the models sampled by the Monte 
Carlo categorized by their "zoology", i.e., whether they 
fit into the category of small-field, large-field, or hybrid. 
We see that all three types of potential are compatible 
with the data, although hybrid class models are preferred 
by the best-fit region. 

Figures El El an d El show the reconstructed poten- 
tials divided by type. Perhaps the only conclusion to 
be drawn here is that the WMAP data places no signif- 
icant constraint on the shape of the inflationary poten- 
tial; many "reasonable" potentials are consistent with 
the data. However, significant portions of the observable 
parameter space are ruled out by WMAP, and future ob- 
servations can be expected to significantly tighten these 
constraints psl l2l| . 

For the particular example of V oc 4 , using Eq. 031 
we find that this choice of potential ruled out to 3<r only 
for N < 40 for the WMAP data set. This constraint 
is even weaker than that claimed by Barger et al. [38|, 
most likely because we allow for a running of the spectral 
index in our constraint. 

When augmenting the WMAP data set with the data 
from seven other CMB experiments, the most noticeable 
improvement in the constraints is a better upper limit to 
the tensor/scalar ratio r, which results in a slightly im- 
proved upper limit on the height of the potential. Also, 
the width of the reconstructed potentials in Planck units 
is somewhat less than in the case of the WMAP-only con- 
straint, showing that the additional data more strongly 
limit the form of the inflationary potential. 
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The combined data rules out a potential with V oc <j) A 
for N < 66 to 3a, effectively killing such models as ob- 
servationally viable candidates for the inflaton potential. 

VII. CONCLUSIONS 

In this paper, we presented an analysis of the WMAP 
data set with an emphasis on parameters relevant for 
distinguishing among the various possible models for in- 
flation. In contrast to previous analyses, we confined 
ourselves to CMB data only and performed a likelihood 
analysis over a grid of models including the full error 
covariance data from the WMAP satellite alone, and in 
conjunction with measurements from BOOMERanG-98, 
MAXIMA-1, DASI, CBI, ACBAR, VSAE, and Archeops. 

We found that the WMAP data alone are consistent 
with a scale-invariant power spectrum, n = 1, with no 
running of the spectral index, dn/dlnk = 0. However, a 
great number of models indicates a compatibility of the 
data with a blue spectral index and a substantial negative 
running. This is consistent with the result published by 
the WMAP team including data from large-scale struc- 
ture measurements and the Lyman-a forest. The WMAP 
result is also consistent with previous CMB experiments. 
The inclusion of previous datasets in the analysis has the 
effect of reducing the error bars and give a better deter- 
mination of the inflationary parameters. Still, no clear 
evidence for the running is present in the combined anal- 
ysis. This result differs from the result obtained in Peiris 
et al. in the case of combined (WMAP+CBI+ ACBAR) 
analysis, where a mild (about 1.5a) evidence for run- 
ning was reported. The different and more conservative 
method of analysis adopted here, and the larger CMB 
dataset used in our paper can explain this difference. 



In addition, we applied the Monte Carlo reconstruc- 
tion technique to generate an ensemble of inflationary 
potentials consistent with observation. Of the three ba- 
sic classes of inflation model, small-field, large-field, and 
hybrid, none are conclusively ruled out, although hybrid 
models are favored by the best fit region. The recon- 
struction process indicates that a wide variety of smooth 
potentials for the inflaton are consistent with the data, 
indicating that the WMAP result is too crude to signif- 
icantly constrain either the height or the shape of the 
inflaton potential. In particular, the lack of evidence for 
tensor fluctuations makes it impossible to constrain the 
energy scale at which inflation takes place. Nonetheless, 
WMAP rules out a large portion of the available param- 
eter space for inflation, itself a significant improvement 
over previous measurements. For the particular case of 
a potential of the form V (</)) — A</> 4 , WMAP rules out 
all such potentials for N < 40 at the 3a level, which 
means that </> 4 potentials are not conclusively ruled out 
by WMAP alone. The combined data set, however, rules 
out 4> 4 models for N < 66, which kills such potentials as 
viable candidates for inflation. 

After this paper first appeared, Leach and Liddle also 
released a reanalysis of the WMAP data |l05j, which is 
in general agreement with our results here. 
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FIG. 1: Regions on the r — n plane. The different types of potentials, small field, large field, and hybrid, occupy different 
regions of the observable parameter space. 




FIG. 3: Models generated by Monte Carlo plotted on the (r, n) plane (black dots). The colored lines are the same models as 
in Fig. 2. For comparison with the models, points are plotted to first order in slow roll. 



14 




FIG. 5: Models generated by Monte Carlo plotted on the (r,dn/d\nk) plane. Points are plotted to second order in slow roll. 
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FIG. 6: Likelihood contours for WMAP alone (left column) and WMAP plus seven other experiments (right column), plotted 
in the (n, dn/dln k) plane (top), the (n, r) plane (center), and the (r, dn/dln k) plane (bottom). The points represent the results 
of the Monte Carlo sampling of inflationary models consistent with la (red, dots) , 2a (blue, triangles) , and 3a (black, stars) 
contours. All points in this and subsequent figures are plotted to second order in slow roll. 
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FIG. 7: Three hundred reconstructed potentials chosen from the sampling shown in Fig. ||3 for the WMAP dataset (left 
column) and the WMAP data set plus seven other experiments (right column). The potentials are color-coded according to 
their likelihoods: la (red) 2a (blue), and 3a (black). The top figure shows the potentials with height and width plotted in 
units of mpi, and the bottom figure shows the same potentials rescaled to all have the same height and width. 
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FIG. 8: Reconstructed potentials from Fig. plotted separately by likelihood: la (red, top), 2a (blue, center), and 3<r (black, 
bottom). The WMAP results are in the left column, and WMAP plus seven other experiments are in the right column. 
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FIG. 10: Likelihood contours plotted in the (n,dn/d\nk) plane (top), the (n,r) plane (center), and the (r,dn/dlnk) plane 
(bottom). The points represent the results of the Monte Carlo sampling, color coded by model type: small-field (red, dots), 
large-field (green, triangles), and hybrid (blue, crosses). The left column is WMAP only while the right column is WMAP plus 
seven other experiments 
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WMAP plus seven other experiments 




FIG. 11: Three hundred reconstructed potentials chosen from the sampling shown in Fig. for WMAP (left column) and 
WMAP plus seven other experiments (right column). The potentials are color-coded according to model type: small-field (red), 
large-field (green), and hybrid (blue). The top figure shows the potentials with height and width plotted in units of m-pi, and 
the bottom figure shows the same potentials rescaled to all have the same height and width. 
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FIG. 12: Reconstructed potentials from Fig. 1111 plotted separately by model type: small-field (red, top), large-field (green, 
center), and hybrid (blue, bottom). The WMAP results are shown in the left column while WMAP plus seven other experiments 
are shown onthe right column. 



